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ABSTRACT 

We report on observations of BL Lacertae during tlie first 1 8 months of Fermi LAT science operations and present results from a 48-day multifrequency 
coordinated campaign from 2008 August 19 to 2008 October 7. The radio to gamma-ray behavior of BL Lac is unveiled during a low activity state thanks to the 
coordinated observations of radio-band (Metsahovi and VLBA), near-IR/optical (Tuorla, Steward, OAGH and MDM) and X-ray (RXTE and Swift ) observatories. 
No variability was resolved in gamma-rays during the campaign, and the brightness level was 15 times lower than the level of the 1997 EGRET outburst. Moderate 
and uncon'elated variability has been detected in UV and X-rays. The X-ray spectrum is found to be concave indicating the transition region between the low and 
high energy component of the spectral energy distribution (SED). VLBA observation detected a synchrotron spectrum self-absorption turnover in the innermost 
part of the radio jet appearing to be elongated and inhomogeneous, and constrained the average magnetic field there to be less than 3 G. Over the following 
months BL Lac appeared variable in gamma-rays, showing flares (in 2009 April and 2010 January). There is no evidence for correlation of the gamma-rays with 
the optical flux monitored from the ground in 18 months. The SED may be described by a single zone or two zone synchrotron self-Compton (SSC) model, 
but a hybrid SSC plus external radiation Compton (ERC) model seems preferred based on the observed variability and the fact that it provides a fit closest to 
equipartition. 

Subject headings: Gamma rays: galaxies - BL Lacertae objects: individual: BL Lac - BL Lacertae objects: 
general - X-rays: galaxies - galaxies: jets - radiation mechanisms: non-thermal 
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1. INTRODUCTION 

BL Lacertae (BL Lac, S4 2200+42, OY 401, B3 2200+420, 
lES 2200+420, 3EG J220 2+4217, IFGL J2 202.8+4216) at 
redshift z = 0.0686 (e.g.. IVermeulen et al.|[T995) was his- 
torically the prototype of the class of active galactic nuclei 
(AGN) for which BL Lac has become the eponym. It is cate- 
gorized as an intermediate synchrotron peaked (ISP) BL Lac 
object. This is based on the latest HSP/ISP/LSP classification 
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proposed by I Abdo e^a P (!201 0a|) to replace th e older HBL- 
IBL-LBL scheme ([Padovani & Giommil 119951: iFossati et al] 
1998). BL Lac occasionally showed peculiar behavior that 
has questioned both this classification and a simple interpre- 
tation of its broad-band emission in terms of synchrotron and 
synchrotron self-Compton (SSC) emission produced by a sin- 
gle blob. In fact, BL Lac is a source that h as shown quite 
compl ex and distinct X-ray spectral behavior ( Madeiski et alj 



1999; Ravasio et alJl2002tlB5ttcher & Reimeiil2004l) and duT 
ing several epochs broad Ha and H/3 emission lines with 
luminosity (^ l O^^ergs"^) com parable to thosejrf type I 
Seyfert galaxies ( IVermeulen et al.lTl995l: [Corbett etan[l996l 
"2000). There is evidence for an increase by ~ 50% in ten 
years of the broad Ha line and an underluminous broad line 
region (BLR), compared to other AGN. The narrow lines and 
radio luminosities of BL Lac match those of of lo w-excitation 
and miniature radio galaxies ('Capet ti et al.ll20 1^. Long term 
radio-optical flux density monitoring and several multiwave- 
length campaigns have been carried out on BL Lac, like the 
past campaigns of th e Whole Earth Blazar Telescope (WEBT, 
IVillata et al.l l2002aHbl) dedicated to thi s source, providing 
very complete datasets (for example. Bloom et al.' '1997t 
^ambruna et al. 1999; Madeiski et al. 1999; Ravasio et aU 
2003; Bottcher et al. 2003; ' Villata et al. 2002a, 2004ali 
12009; Raiteri et al.. 2009, 2010.; Marscher et al.„2008) . In par- 
ticular in summer 1997 BL Lac showed the largest optical 
outburst ev er recorded in almost 33 years dNesci et alJI 19981; 
iTosti et alj ri999; Villata et al. 2004b). 

In the X-ray band the two broad components of the spectral 
energy distribution (SED) are overlapping and the radiation 
at this band is at times dominated by the high-energy end of 
the synchrotron emission, while at other occasions it is dom- 
inated by the low-frequency portion of the high-energy com- 
ponent (e.g., 'Maddslaitanillll Eva^^ The 
rapid X-ray variability is mainly restricted to the low-energy 
excess portion of the X-ra y spectrum, presumably pr oduced 
by synchrotron radiation dRavasio et alj 120021 l2003h . The 
X-ray variability of BL Lac, except fo r a few major flaring 
events, has a log-normal distribution (Gieb els & Degrarigg 
12009), meaning that the emission is a multiplicative product 
of a large number of independent random events. 

Gamma-ray observations by EGRET resulted in several de- 
tections (Figure [Til. Observations after 1995 resulted in an 
average 7-ray flux above 100 MeV of (40 ± 12) x 10~^ 
ph cm"2 s"i (ICatanese etalJll997h and in a flare in 1997 
at a level (171 ± 42) x 10"^ ph cm'^ with 10.2cr sig- 
nificance (Bloom et al] 119971) . Correlated 7-rays and opti- 
cal flaring emission were observed during the EGRET era 
(iBloom et alJll99'7h . Very high energy (VHE) 7-ray emis- 
sion was claimed by the Crimean GT-48 (iNeshpor et al.l20()l]) 

School of Pure and Applied Natural Sciences, University of Kalmar, 
SE-391 82 Kalmar, Sweden 

" Tuorla Observatory, University of Turku, FI-21500 Piikkio, Finland 
Department of Physics and Astronomy, Ohio University, Athens, OH 
45701, USA 

Instituto Nacional de Astrofisica, Optica y Electronica, Tonantzintla, 
Puebla 72840, Mexico 

Instituto de Astronomia y Meteorologia, CUCEI, Universidad de 
Guadalajai'a, 44130 Guadalajara , Jalisco, Mexico 

Aalto University Metsahovi Radio Observatory, FIN-02540 Kylmala, 
Finland 

^ Institute for Astrophysical Research, Boston University, Boston, MA 
02215, USA 

Finnish Centre for Astronomy with ESO (FINCA), University of 
Turku, FI-21500 Piikkio, Finland 



The first Fermi multifrequency campaign on BL Lac 



3 



< — 1 




[ 1 1 










: EGRET 






> ■ 


E 


150 


L (E>100MeV) 






o 




















CL 


100 








CO 












O 






1 




X 


50 








k : 







■,tI,,„„x.l.t....T 









1 99 1 1 992 1 993 1 994 1 995 1 996 1 997 1 998 



Fig. 1. — Gamma-ray light curve of BL Lac obtained by EGRET in the 
period 1991-1997. Data points have source detection significance above 4a 
iNandikotkur et aU2007h . 

and HEGRA (iKranichI 120031) atmospheric Cherenkov tele- 
scopes. A significant (> 5 cr) detection above 200 G eV was 
made only by the MAGIC telescope in 2005 ( Albe rt et all 
120071) . Gamma-ray emission was explained with the require- 
ment of Comptonization of external-jet photons (external- 
radiat ion Compton, ERC) in addition to the SSC (in- j et) emis- 
sion (MidiikiitiDIIlS IIMcEE^^ Super- 
luminal motion of /3aa D up to (10 . 57 zt 0.74) has been ob - 
served in this object (Lister etal] 120091: iDenn et alj |2000). 
Following results of past multifrequency campaigns, possibly 
distinct VLBI jet structures are assumed to contribute, some- 
times with delays, to the radio flux density light curves, and 
sometimes are suggested to be responsible for optical to TeV 
7-ray flares (for example, iBach et aH l2006t iMarscher et alJ 
l2008h . 

In this paper, we report on LAT observations of BL Lac 
during the first 1 8 months of Fermi science operations (from 
2008 August 04 to 2010 February 04). The LAT's reasonably 
uniform exposure, high sensitivity, and continuous sky moni- 
toring make it an excellent instrument around which to orga- 
nize a simultaneous multifrequency campaign. The so-called 
Fermi planned intensive campaign (PIC) dedicated to BL Lac 
was performed in 2008 August 19 - 2008 October 7 (MJD 
54697.8 - 54746), during roughly the first two months of sci- 
ence operations, irrespective of the brightness of the source. 
This was part of a series of Ferm/-LAT-collaboration PICs 
(ffosti 2007; Thompson 2007), involving observing proposals 
submitted to the RXTE (Cycle 12) and Swift (Cycle 4) X- 
ray satellites, and organized in advance of the Fermi launch. 
This ensured access to the facilities allowing the best multi- 
wavelength (MW) coverage. The aim of the campaign is to 
shed light for the first time on the broad-band radio-to-7-ray 
SED, including the high energy (X-ray and 7-ray) behavior, 
during a low activity phase of the source. The 18-month LAT 
light curve shows that BL Lac was variable in 7-rays for most 
of these 18 months, with the exception of the first 2 months, 
which corresponds to the period of the MW campaign. In sec- 
tion|2]the light curve by Fermi -LAT during the first 18 months 
of survey and simultaneous optical long-term monitoring data 
are presented joined to the analysis of the 7-ray spectrum. In 
section|3]the multiwaveband radio-band to X-ray flux density 
and parsec-scale radio structure observations collected during 
the campaign period are reported and discussed. The assem- 
bled spectral energy distribution and modeling are reported in 



section Upl and conclusions in section |4] 

A ACDM cos mology with values given within Icr of the 
WMAP results jKomatsu et al.1 1200 9*), namely flm = 0.27, 
r^A = 0.73, and Hq = 71 km s"^ Mpc^^ is used. 

2. GAMMA-RAY OBSERVATIONS AND RESULTS BY FERMl-LAT 
2.1. LAT analysis and observations 

The Large Area Telescope (L AT), on b oard the 
Fermi Gamma-ray Space Telescope ( I Atwood ei t al. 2009), 
is a pair-conversion 7-ray telescope, sensitive to photon 
energies from about 20 MeV up to >300 GeV. It consists 
of a tracker (composed of two sections, front and back, 
with different angular resolutions), a calorimeter and an anti- 
coincidence system to reject the charged-particle background. 
The LAT has a large peak effective area (r-^ 8000 cm^ for 
1 GeV photons in the event class considered here), viewing 
« 2.4 sr of the sky with angular resolution (68% containment 
radius) better than « 1° at £^ — 1 GeV. The lai-ge field of 
view, improved effective area and sensitivity and the survey 
nominal mode make Fermi-LAJi an optimal all-sky hunter 
for high-energy flares and an unprecedented monitor of 7-ray 
sources. The data set used in this paper was collected during 
the first 18 months of Fermi science obserations, from 2008 
August 4 to 2010 February 4 (about 550 days, 78 weeks from 
MJD 54682.7 to 55232.9, as shown in the weekly fight curve 
in Figure |2] upper panel). This interval includes the period 
chosen for the first Fermi MW campaign on BL Lac (2008, 
Aug. 19 - Oct. 7, MJD 54697.8-54746, i.e. about 48 days, 
indicated by the line in Figure |2] upper panel). Fermi -LAT 
data analysis was performed with the standard Fermi -LAT 
ScienceTools software package^ using version v9r 1 5p5 . 
Only events belonging to the "Diffuse" class in the energy 
range 0.1 — 100 GeV were used in the analysis. Instrument 
response functions (IRFs) used were P6_V3_DIFFUSE. In 
order to provide protection against significant background 
contamination by Earth-limb 7 rays, all events with zenith 
angles > 105° were excluded. 

The 18-month light curve was built using 1-week time 
bins (Figures |2] and |4|i. For each time bin the inte- 
grated flux (E>100 MeV) values were computed using the 
maximum-likelihood algorithm implemented in the science 
tool gtlike. For each time bin we analyzed a region of 
interest (Rol) of 12° in radius, centered on the position of the 
source. All point s ources listed in the IFGL (1-year) LAT 
catalog (lAbdo et al.ll201 0b) within 19° from BL Lacertae and 
having test statistic^ TS > 50 and fluxes above 10^^ ph 
cm^^ s~^ were included in the Rol model using a power- 
law spectrum (dN/dE oc E~^, where F is the photon index). 
The isotropic background (the sum of the residual instrumen- 
tal background and extragalactic diffuse 7-ray background) 
was included in the Rol model using the standard model 
file isotropic_iem_v02 . txtPl and the Galactic diffuse 
emission was included in the modeling using the file and the 
standard file gll_iem_v02 . fit. In the final light curve 
computation the photon index value was frozen to the value 
resulting from the likelihood analysis on the entire period. For 
each time bin, if the TS value for the source was TS < 4 or 
the number of model predicted photons Npred < 10, the value 

^^ |http : //f ermi . gsf c ■ nasa ■ gov / ssc/data/analy sis/documentation 
The test statistic is defined as TS = 2Alog(Likelihood) between mod- 
els with and withou t the source and it is a measure of the source significance 
)MattoxetaU199eh . 

http ; //f ermi . gsf c .nasa ■ gov/ ssc/data/access/lat /BackgroundM 
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Fig. 2. — Top panel: 18-month weekly integrated flux (E > lOOMeV) light curve of BL Lac obtained obtained from 2008, Aug. 04 to 2010, Feb. 04. 
Bottom panel: superposed optical V-band and R-band observed flux (not corrected for interstellar absorption and host galaxy contribution) light curves obtained 
respectively at the Steward Observatory, Kitt Peak, Arizona, USA, and Tuoria Obs 



of the fluxes were replaced by the 2-a upper Hmits. AH errors 
reported in the figures or quoted in the text are l-cr statistical 
errors. The esti mated relative syste matic uncertainty on the 
flux, according to lAbdo et alj (12010b.) and reflecting the rela- 
tive systematic uncertainty on effective area, is set to 10% at 
100 MeV, 5% at 500 MeV and 20% at 10 GeV. 

2.2. Gamma-ray spectral and temporal behavior 

The 7-ray spectral analysis of BL Lacertae was performed 
both for the entire 18-month period, and for the 48 days of 
the campaign. These spectra are shown in Figure [3] The flux 
value for each energy bin and for the two periods are reported 
in Table [T] The most energetic photon observed within the 
95% point spread function containment radius of the source 
for the 18 -month period had energy 70 GeV, while for the 
campaign period the highest energy photon had 10 GeV. The 
energy-spectrum binning was built requiring TS > 50 and/or 
model-predicted source photons > 8, except for the last 
bin, which had an upper limit rather than a detection. The 
18-month spectrum results in an integrated average flux of 
(19.85 ± 0.40) X IQ-** ph cm^^ g-i , in the 0.1 - 100 GeV 
range, with TS ^ 2375 and photon index T = 2.38 ± 0.01, 
while the value for the 48 -day interval of the multifrequency 
campaign, (11.5 ± 2.7) x 10^^ ph cm^^ s^^ , and the photon 
index slightly flatter (T = 2.27 ± 0.10), with nTS ^ 111. 

In order to quantify the departure of the 18 -month spectrum 
from a power law shape, we use a spectral curvature index C 



□ LAT spectrum campaign period 
—Power-law 

• LAT spectrum 18 montiis 
' LogParabola 




lO-* 10" 
Energy [MeV] 

Fig. 3. — 18-month averaged energy spectrum of BL Lac obtained with 
a band likelihood over the accumulated observations from 2008 Aug. 04 to 
2010 Feb. 04 (MJD 54682.7-55232.9). The spectrum shows a departure from 
a simple power law that can be fitted by a log parabola function. The blue 
(open square) symbols represent the spectrum corresponding only to the 48- 
day period of the multifrequency campaign. 
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TABLE 1 

BL Lac 7-ray spectrum obtained with likelihood 

ANALYSIS IN EACH ENERGY BAND. 





7-ray flux 


Syst. 


[McV] 


[xlO"**ph cm^^s"^] 




1 8 months 


100-200 


7.62±1.03 


10% 


200-562 


5.14±0.27 


5% 


562-794 


0.73±0.06 


5% 


794-1122 


0.45±0.04 


5% 


1122-1585 


0.31 ±0.03 


10% 


1585-2239 


0.18±0.02 


10% 


2239-3162 


0.11 ±0.02 


10% 


3162-6310 


0.10±0.01 


15% 


6310-35481 


0.04±0.01 


15% 


Aug 19 - Sep 9 


100-562 


7.18±2.12 


5% 


562-2239 


1.28±0.24 


10% 


2239-12589 


0.17±0.07 


15% 



(4) 



where i is the index for a particular energy bin, Fi is the ob- 
served flux, Ff^ is the flux in the ith bin predicted by the 
global power-law fit, a is the error in Fi, and /['^^ is the rel- 
ative systematic uncertainty in that bin (reported in Table [T]i. 
The power-law fit has two parameters: the normalization and 
photon index. The curvature parameter is expected to follow 
a distribution with 9 — 2 = 7 degrees of freedom (d.o.f.) 
if the power law hypothesis is true. With 99% confidence 
the spectral shape is significantly different from a power law 
when C > 18.48. For the accumulated 18-month data on BL 
Lac we found C = 20.88 (0.1 < E < 100 GeV), which im- 
plies a curved spectrum. The same procedure, when applied 
to the spectrum of the 48-day multifrequency campaign, re- 
sulted in C = 2.83 (0.1 < E < 100 GeV), implying there is 
no significant evidence for deviation from a simple power-law 
in this spectrum (Fi gures [3] and [TOb . This does not necessar- 
ily mean that there is no curvature during this 48-day period, 
because it could be due to the lower emission state and accu- 
mulated statistics. 

To further explore curvature, the 18-month spectrum was 
fit with a log-parabola function in the 0.1 — 300 GeV energy 
band. This was performed in two runs: first a fit with Ehreak 
left free, then another with Ebreak frozen at the value found 
in the first run (300 MeV) to improve the calculation of the 
other parameter values. The result of the latter fit gives 

dN/dE={2m ± 0.07) X 10"^° 

X (£;/300 MeV")"(^^-^^=^"-°^^+(°-"^^"-°^^'°'5(£;/300 MeV)) 

with TS = 2384.1. The higher TS value shows that the 18 
months LAT spectrum does not follow a simple power law 
function, and the log-parabola model can describe this spec- 
trum. 

The weekly light curve in the 0.1-100 GeV band for the 18- 
month period is presented in Figure |2] (upper panel). In order 
to quantify the vari ability, the variability index as defined in 
lAbdo et al.l(l20T0bh . was computed from 



V = Y,w,{F,-F^ 



where 



1 



Wi 



„2 + [Jrelp^Y 



(2) 



(3) 



and the index i runs over all data points except the upper limits 
(61 bins). Here Fi and cji are flux values and the relative 
statistical errors, respectively, and f^-ei is equal to 3% of th e 
flux for each interval, as sugg ested by lAbdo et all ( l2010bl) . 
For the 18 months, F^t was found to be 17.48 x 10~* ph 
cm~^ s^^ and V = 193 with 60 d.o.f., which is consistent 
with variability (the probability the source is non-variable is 
less than 10~^^). The same values computed just for the first 
9 weeks, containing the period of the MW campaign, give 
F^t = 10.67 X 10"* ph cm-2 g-i and V = 0.6 with 8 
d.o.f., which is consistent with a non-variability hypothesis at 
> 99% confidence. Each of the LAT flux data points in the 
upper panel of Figure H has TS > 10. 

This interval with no observed 7-ray variability character- 
izes the intensive coordinated campaign period. After the 
campaign, starting from about halfway through 2008 Novem- 
ber until 2010 February 4, BL Lac showed variable weekly 
7-ray flux with V = 162 with 51 d.o.f., F^t = 19.15 x 10"^ 
ph cm^^ s^^ . In particular, a 1-week 7-ray flare occurred in 
the week MJD 54928.5-54935.5 (2009 Apr 7-14), with a flux 
value of 72 x 10^^ ph cm^^ s^^ . This is about 4 times greater 
than the mean flux value (F^t) computed above. Moreover 
this is the highest weekly flux yet detected from BL Lac, al- 
though another likely flare was dete cted in th e week around 
MJD 55221 (2010 Jan. 25. iSokolovsky et al]l2010al) . After 
the 2009 April flare the source also showed slightly higher 
activity and flux. 

The 78 points of the 18-month hght curve were also used 
to com pute the globa l normalized excess variance, as de- 
fined inMuihanltaD(l2Q0l;|Mdoetd^ a%xs = 
(5*^ — a1„)/x^, where S'^ is the variance of the light curve 
and dg^^ = erf + crlys, the sum in quadrature of the statisti- 
cal uncertainty of the flux in the time bin and the systematic 
error estimate (0.03(Fi)). This quantity measured an intrinsic 
variability amplitu de of 0.23 ib 0.03, in good agreement with 
the value found in lAbdo et aP (l2010dh reporting the weekly 
light curve of BL Lac (flux E > 300 MeV) during the first 1 1 
months. This confirms again that BL Lac was variable after 
the first couple of months of Fermi science operations. 

3. OBSERVATIONS AND RESULTS FROM THE MULTIFREQUENCY 
CAMPAIGN 

The 48-day PIC on BL Lac involved the participation of 
the RXTE and Swift X-ray satellites, and of ground-based ra- 
dio and optical observatories. These included the Metsahovi 
13.7 m radio-telescope operating at 37 GHz in Finland; the 
Very Long Baseline Array (VLBA) in the USA, which took a 
multi-waveband flux-structure observation on 2008 Septem- 
ber 2; two telescopes of the Tuorla Observatory, Finland, and 
two telescopes of the Steward Observatory, USA, for the long- 
term and single-band optical monitoring as presented in Fig- 
ure|2l and the 2.1 m optical telescope of the Observatorio As- 
troffsico Guillermo Haro (OAGH, Mexico) operating in near- 
infrared, and the 1.3 m McGraw-Hill optical Telescope of the 
MDM Observatory (Arizona, USA), for further multi-band 
optical snapshots during the campaign period (Table |2]l. 

3.1. X-ray and UV observations and results 

The RXTE/PCA space observatory performed sub-daily 
monitoring of BL Lac for 20 days, during the same epoch 
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Fig. 4. — Simultaneous multifrequency flux light curves of BL Lac from the 48-days Fermi intensive campaign obtained with Fermi -LAT. Swift - XRT and 
UVOT, OAGH, and Metsahovi observations. LAT weekly fluxes reported in the upper panel for this period are all TS > 10 detections. UV-to-near-IR fluxes are 
corrected for absorption. 



as LAT and Swift, for a total of 80 pointings between 2008 
August 20 and September 9. Only PCU2 was active, for a to- 
tal exposure of 157.3 ks. The PCA STANDARD2 data were 
reduced and analyzed with the routines in HEASOFT V6.8 
using the filtering criteria recommended by the RXTE Guest 
Observer Facility. Only the top-layer events were processed, 



and a check was performed in the 40-100 keV range to assure 
that the model background reproduced the observed one. The 
average net count rate in the 3-18 keV band is 0.55 ± 0.01 
ctss^^pcu"^. The RXTE spectra were extracted and fitted 
separately for each pointing, and summed to obtain the aver- 
age spectrum. Each spectrum is well-fit by a single power- 
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TABLE 2 

BL Lac coordinated multifrequency campaign (PIC) 

AND LONG-TERM MONITORING OBSERVATIONS 



PIC (48 days) 








Instrument 


Energy Range 


2008 Epoch range 


# obs. 


VLBA 


4.6 - 43.2 GHz 


Sep 2 


7 


Metsiihovi 


37 GHz 


Aug 20 - Oct 6 


22 


OAGH 


JHK 


Sep 6 - Oct 6 


18 


MDM 


UB VRI 


Oct 6 - 10 


15 


Swift -UVOT 


W2 M2 WI U B V 


Aug 20 - Oct 2 


141 


Swift -XRT 


0.4-8 KeV 


Aug 20 -Sep 18 


24 


RXTE -PCA 


3-18 KeV 


Aug 20 - Sep 8 


19 


Fermi-LAT 


100 MeV - 100 GeV 


Aug 19 - Oct 7 


48 days 


1 8 months 








Instrument 


Energy Range 


Epoch range (MJD) 


#obs. 


Tuorla 


R 


54709.8 -5519L8 


162 


Steward 


V 


54743.2-55213.1 


89 


Fermi-LAT 


100 MeV - 100 GeV 


54682.7 - 55070 


78 weeks 



law model with Galactic absorption. Th e Galactic Nn wa s 
fixed at two values: that measured by lElvis et aTl (119891) . 
based on dedicated 21-cm observations (2.015 x 10^^ cm~^); 
and the sum of this value wi th that inferred from millimeter 
observations (iLucas & Lisztj[T993 ). which includes the con- 
tribution from molecular hydrogen along the line of sight 
(Nn ~ 3.6 X 10^^ cm-2). In the RXTE band, however, 
this difference yields negligible effects, with a difference in 
spectral indices ATx ^ 0.02. 

The RXTE data showed modest flux variations, with a rate 
which remained constant at ^ 0.6 cts s^^ up to MJD 54708, 
and slowly decreasing to ^ 0.4 cts s~^ towards MJD 54718. 
Spectral variations were modest as well, ranging from 1.95 ± 
0.09 to 2.28 ± 0.07. The total average spectrum is well-fit 
with a single power-law model, with Fx — 2.18 ± 0.06, and 
an unabsorbed flux in the 2-10 keV band of F2-1Q — (7.2 ± 
0.3) X 10-12 ergcm-^s-i and xl = 0.81 (15 d.o.f.). This 
spectrum can also be fit by a concave log parabola or broken 
power-law model, with the slope hardening above 8-9 keV 
from ^ 2.2 to ~ 1.9 (xf. — 0.7); however, the low statistics do 
not allow one to determine if the improvement is significant, 
and the single, rather flat (Fx — 2, therefore i/Fi, oc ly'^), 
power law gives a good fit. 

In the context of the overall SED, the low X-ray flux and 
the Fx — 2 spectrum (flat X-ray SED) is evidence that the 
RXTE band might correspond to the passage between the two 
humps of the blazars' SED, as typical for ISPs. In the fol- 
lowing we show that the Swift -XRT results provide us further 
evidence for this conclusion. 

The Swift gamma-ray burst satelhte (iGehrels et al.ll2004h 
performed a daily monitoring simultaneous to RXTE and LAT 
observations for 20 days (from 2008 Aug. 20, 15:19 UT, to 
2008 Sep. 09, 02:37 UT; MJD 54698.64 - 54718.11), plus 3 
more separated observations after 18 September. Data reduc- 
tion and analysis was performed running a processing script 
customized for the XRT and UVOT data (Donato et al., in 
preparation). The script reprocesses the Swift data, stored in 
the HEASARC archive using the standard HEASoft software 
(version 6.8) and the latest calibration database (20091 130 for 
XRT and 20100129 for UVOT). The reduction of XRT data 
consists of running xrtpipeline, selecting only the events with 
0-12 grades in photon counting mode (PC). The UVOT im- 
age mode data are processed following the steps reported in 
the UVOT Software Guide 2.2. 
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Fig. 5 . — Simultaneous near-IR to IJV SEDs of BL Lac at different epochs, 
collected with UVOT ( Aug. 20, 23, Sept. I, 4), OAGH (Sept. 6), and MDM 
(Oct. 6) observations. 



The XRT data analysis is performed using the script 
xrtgrblc, available in the HEASoft software package. In 
brief, the script selects the best source and background ex- 
traction regions based on the source intensity (a circle and an 
annulus, respectively). In the case of BL Lac, the source ex- 
traction region has typically a 55"radius and the background 
region has a 1 10"-210"inner-outer radius. Field sources are 
excluded applying circular masks whose radius depends on 
the field source intensity. Using these regions, source and 
background count rates, spectra and event lists are extracted. 
The net source count rate must be corrected for irregularities 
in the exposure map and the Point-Spread Function (PSF). 
The total correction factor is obtained using the HEASoft tool 
xrtlccorr. 

The script that handles the UVOT analysis is uvotgrblc 
(available within HEASoft). It determines the presence of 
field sources and excludes them from the background region 
with circular regions whose sizes depend on their intensity. 
The optimal background region is chosen after comparing 3 
annular regions centered on the main source in the summed 
images. The region with the lowest background is selected. 
The inner and outer radii are 27"and 35", respectively, for 
all the filters except V, for which the two radii are 35"and 
42". The source extraction region is also intensity dependent: 
The script selects the aperture size between 3"and 5"based 
on the observed source count-rate and the presence of close 
field sources. BL Lac is a relatively bright object and an 
extraction region of 5"is preferred. The field of view in 
the V and B filters is full of field sources and in particular 
there is a bright star east of the blazar To avoid contami- 
nation from such a source, the extraction region is reduced 
to 3". The script estimates the best source position using 
the task uvotcentroid and the photometry using the task 
uvotsource. The obtained values are aheady corrected for 
aperture effects while we used the dust maps of Schlegel et al] 
(1998) and the Milky Way extinction curve of Pei (1992) to 
compensate the Galactic extinction. XRT and UVOT light 
curves and SEDs are shown in Figures |4] |5] and|6] 

As with the harder energy band with RXTE, BL Lac does 
not show any major sign of variability in the Swift-XRT data. 
The same Galactic (TVh — 3.6 x 10^^ cm~-^) fixed value 
was used also for the XRT analysis. The 2-10 keV flux de- 



8 



Abdo et al. 



BLLacertae, RXTE total 

PCU2, power-law 



1 1 1 1 1 1 1 








++. , , 


++ 






1 





5 10 20 

Energy (keV) 



Swift XRT 0.3-1 OkeV 




0.5 1 2 5 

Energy (keV) 



Fig. 6. — Top panel: RXTE combined 3-18 keV spectrum of BL Lac ex- 
tracted cumulating the 20 days of observations of this X-ray satellite during 
the campaign. Bottom panel: Swift -XRT combined 0.4-8 keV spectra of 
BL Lac extracted cumulating the 24 days of observations of this X-ray satel- 
lite during the campaign. 

rived from Swift observations lies between 6 and 8 xlO^^^ 
erg cm^^ s^^with a hint of a lower intensity at the end 
of the campaign when the flux was below 5 x 10^^^ erg 
cm^^ s^^. We found that in the day-by-day spectral anal- 
ysis the photon index does not change significantly either. 
Using a single power-law model with fixed Galactic absorp- 
tion, Leiden/ Argentine/Bonn (LAB) Survey of Galactic HI 
dKalberla et al. 2005) weighted average equal to 1.73 x 10 
cm~^), the slope varies between 1.75 and 2.06. Since the typi- 
cal error on this parameter is of the order of 0. 1 , the variability 
is within 2 sigma. 

The absence of spectral changes allows us to extract an ac- 
cumulated spectrum. As a first step, a comparison of the XRT 
and RXTE spectra accumulated over the same period and us- 
ing the same lower energy threshold of 3 keV results in very 
good agreement with the power-law photon index Tx found 
by RXTE . The best fit in the 3-8 keV range is obtained with a 
power-law whose slope is 2.19lo;i[5. The 2-10 keV integrated 



flux is 6.37to\l X 10"^^ erg cm^^ g-i x-ray photon 
index is consistent within la with the results obtained using 
RXTE data. 

Using the entire XRT energy band from 0.3 to 10 keV (Fig- 
ure |6]l and profiting from the increased statistics (the spec- 
trum is binned with 100 counts/bin), the spectral fit is more 
constrained and underlies two emission components: the best 
fit (xl — 0.82 with 146 d.o.f.) is obtained with a broken 
power-law with the energy break located at 1.71^q']^^ keV 
and fixed absorption at the Galactic value. The soft and hard 
photon indices are 

2.57t°:o5 and 2.09±0.05, respectively and 
therefore the break in this XRT spectrum is significant com- 
pared to the power law. The 0.3-10 keV integrated flux is 
9-67to.io X 10-12 ergcm-2 g-i. 

The UVOT optical-UV light curves showed some variabil- 
ity during the campaign period (FigurelH, but this is not well 
correlated with the X-ray or 7-ray emission (which is not vari- 
able). The daily simultaneous near-IR, optical, and UV SEDs, 
reconstructed thanks to UVOT, and OAGH observations (re- 
ported in Figure lU, showed no significant spectral variability 
and a trend consistent with a single power law except for the 
2 higher frequency UV bands. This may be in agreement with 
the hypothesis of UV excess due to therm al emission from an 
accretion disk jRaiteri et al.ll2009ll2010l) . 

3.2. Radio band flux-structure observations and results 

As part of an ongoing blazar monitoring program, the 
37 GHz observations were made with the 13.7 m diame- 
ter Metsahovi radio telescope, which is a radome enclosed 
paraboloid antenna situated in Finland. The measurements 
were made with a 1 GHz-band dual beam receiver centered at 
36.8 GHz. The high electron mobility pseudomorphic transis- 
tor front end operates at room temperature. The observations 
are ON-ON, alternating the source and the sky in each feed 
horn. A typical integration time to obtain one flux density data 
point is 1200-1400 s. The detection limit of the telescope at 
37 GHz is on the order of 0.2 Jy under optimal conditions. 
Data points with a signal-to-noise ratio < 4 are handled as 
non-detections. The flux density scale is set by observations 
of DR 21. Sources 3C 84 and 3C 274 are used as secondary 
calibrators. A detai led description o n the d ata reduction and 
analysis is given in iTerasranta et al.l (I1998I) . The error esti- 
mate in the flux density includes the contribution from the 
measurement rms and the uncertainty of the absolute calibra- 
tion. 

BL Lac was also observed with the Very Long Baseline Ar- 
ray (VLBA) at seven frequencies (4.6, 5.0, 8.1, 8.4, 15.4, 23.8 
and 43.2 GHz) in the framework of a surv ey of parsec-scale 
radio s pectra of twenty 7-ray bright blazars (ISokolovskv et alj 
l2010bl) . The multifrequency VLBA observation was carried 
out on 2008 September 2. The data reduction was conducted 
in the standard manner using the AIPS package (iGreisenI 
Il990h . The final amplitude calibration accuracy is estimated 
to be ~ 5 % at 4.6-15.4 GHz rang e and ^ 10 % at 23. 8 and 
43.2 GHz. The Difmap software ( Shepherd et"anil994j) was 
used for imaging and modeling of the visibility data. 

The VLBA images (Fig. |7] and [8]) reveal a wide, rather 
smooth, curved jet extending ^ 50 mas (at 5 GHz) south-east 
from the bright compact core. A few distinct, bright emis- 
sion features aligned along the south-southwest direction can 
be seen in the inner jet at higher frequencies. We have mod- 
elled the observed brightness distribution with a small num- 
ber of model components having two-dimensional Gaussian 
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Fig. 7. — Inner jet of BL Lacertae as observed by the VLBA on &pteim)er 02, ^Ti()f5. Irnages'at different frequencies are shifted by 8 mas in relative Right 
Ascension. For the 15.4, 23.8 and 43.2 GHz image map peaks are 1.69, 1.52, 1.32 Jy/beam, first contours are 1.70, 3.00, 5.00 mJy/beam respectively. The contour 
levels are increased by a factor of 3. Beam size (natural weighting) for each frequency is indicated by the cross to the left of the con'esponding image. 



■2 I 



2200-1-420 



2008-09-02 




\ 1 



i 



30 



20 10 -10 

Relative Right Ascension (milliarcsec) 



-20 



Fig. 8. — Spectral index map (Fi, ~ a is shown in color) of BL Lac- 

ertae constructed using VLBA observations at 4.6, 5.0, 8.1 and 8.4 GHz. The 
overlaid contours represent total intensity at 8.4 GHz (the peak intensity is 
1.75 Jy/beam, the first contour is 0.70 mJy/beam, the beam size is 1.57 X 1.22 
mas at PA 10.4°.). The spectral index map was smoothed by a median filter 
with radius equal to the indicated beam size. The map shows an optically thin 
jet with a ~ 0.7 and the self-absorbed core region (a < 0.0). A 2D cross- 
con'elation technique using the optically thin part of the jet was employed to 
align images at different frequencies allowing reliable extraction of the spec- 
tral information. The spectral steepening towards the jet edges visible on the 
spectral index map occurs on the angular scale comparable to the beam size 
and is likely an effect of the uneven MD-coverage at different frequencies. 



profiles. The core region (i.e., the bright feature at the north 
end of the jet) is elongated in the 43 GHz image and can be 
modelled by two Gaussian components separated by 0.25 mas 
(0.32 pc projected distance). At a distances of 1.5 mas and 3.4 
mas from the core, there are two other distinct emission fea- 
tures in the jet. A table listing the paramet ers of the distinct 
features ("components") can be found in ISokolovskv et alj 
(I2010bh . 

The inner 0.25 mas part of the jet (i.e., the "core") may 



contain two (or more) distinct emission features or it may be 
a continuous emission region. The angular resolution even 
at 43 GHz is not sufficient to distinguish between these pos- 
sibilities. However, it is evident that the radio spectrum is 
changing along this region and it cannot be described by a 
single, uniform, self-absorbed, synchrotron emitting compo- 
nent. A turnover caused by the synchrotron self-absorption 
is detected in the averaged core spectrum at a frequency of 
^ 12 GHz (an average over the whole area with the nega- 
tive spectral index at Fig. [U. Due to the inhomogeneity, the 
innermost component at mm-wavelengths most likely has an 
even higher t urnover frequency, > 40 G Hz. Using the method 
described in lSokolovskv et al.l ( l2010bl) . the magnetic field, B, 
of the core, in the frame of the relativistic jet, can be con- 
strained, given the Doppler factor, T> = (r[l — /3 cos9])~^ , 
where F = (1 — and 9 here are the bulk Lorentz 

factor and jet angle to the line o f sight, respectively . Assum- 
ing a Doppler factor V = 7.3 (iHovatta et al.1 1200 9). an up- 
per limit can be placed on the magnetic field strength in the 
core: B < 3G. We note that this upper limit corresponds to 
a typical value across an extended and inhomogeneous region 
contributing the bulk of the radio emission in 4.6 - 43 GHz 
range. The magnetic field strength can exceed the above limit 
locally or in a more compact regions hidden from sight in the 
observed frequency range by synchrotron opacity. 

3.3. Optical and near IR observations and results 

During the campaign several multi-band near-infrared and 
optical observations of BL Lac were also performed from 
the ground by the 2.1 m telescope of the OAGH observa- 
tory, Sonora, Mexico, and the 1.3 m McGraw-Hill Telescope 
of the MDM Observatory, Arizona, USA. In addition a long- 
term monitoring in optical R and V bands during the first 18 
months of Fermi observations (Figure |2]l, was performed by 
the observing monitoring programs at the Tuorla Observatory 
(1 m Tuorla and 0.35 m KVA telescopes), Turku, Finland, and 
at the Steward Observatory (2.3 m Bok and 1 .5 m Kuiper tele- 
scopes), USA respectively. 

Optical (UBVRI) data were taken with the 1.3 m 
McGraw-Hill Telescope (MDM observatory). Exposure times 
ranged from 40 s (R band) to 120 s {U band). The raw data 
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were bias and flat-field corrected using IRAF. Instrumental 
m agnitudes of BL La c and comparison stars B, C, H, and K 
of lSmithetal] ( 119851) were extracted using DAOPHOT and 
subsequently converted to physical magnitudes through dif- 
ferential photometry. 

INAOE (Instituto Nacional de Astrofisica, Optica y 
Electronica) operates the Observatorio Astroffsico Guillermo 
Haro (OAGH) located in the Mexican state of Sonora. The 
2.1m telescope of the OAGH has a current allocation of over 
50 nights per semester dedicated to the study of 7-ray sources, 
through optical photometry and spectroscopy and infrared 
photometry. The Cananea Near Infrared Camera, CANICA, 
is equipped with a Rockwell 1024 x 1024 pixel Hawaii in- 
frared detector working at 75 K and standard J, H and 
bands. The plate scale in CANICA is 0.32 arcsec per pixel. 
Observations are usually carried out in series of 10 dithered 
frames in each filter. Data are co-added after bias and flat- 
field corrections. CANICA observations of BL Lac presented 
in this work were made on seven epochs between MJD 54715 
and 54747. The data are shown in Figures|4]and|5] and do not 
show signs of variability. 

The photometric and polarimetry monitoring program at the 
Tuorla Observatory, a division of the Department of Physics 
and Astronomy at the University of Turku, Finland, is per- 
formed through the 1 .03 m telescope at Tuorla Observatory, 
and the 35cm telescope at the KVA observatory on La Palma, 
Canary islands, Spain. Most of the monitored sources are BL 
Lac objects listed in Costamante & Ghisellini (2002), as po- 
tential TeV 7-ray sourceS All phot ometric measures a re 
taken using the R-filter (other details in iTakalo et al.ll2008l) . 

The optical spectropolarimetry and spectrophotometry of 
blazars at the University of Arizona, is performed on a regular 
basis since the launch of Fermi using the Steward Observa- 
tory 2.3 m Bok telescope on Kitt Peak, Arizona, and 1.54 m 
Kuiper telescope on Mt. Bigelow, Arizona, USA. Depend- 
ing on weather conditions, BL Lac is observed nightly dur- 
ing each of the monthly monitoring campaigns, which typi- 
cally last a week. All of the measurements have been made 
with the SPOL specti-opolarimeter ( Schmidt et alJll99Z) . The 
Steward Observatory blazar monitorin g program associated 
with the Fermi mission is described bv lSchmidt et al.l (120091) 
and the data are publicly availabl4Z3- Details concerning the 
reduction and calibration of both the polarime try and photom- 
etry can also be found in Schmidt et al.l (12009) and references 
therein. The data products include high signal-to-noise flux 
and linear polarization spectra spanning A = 4000 — 7550A 
and differential V-band photometry. 

3.4. Multifrequency correlation and SED 

The optical R- and V^-band light curve in the bottom panel 
of Figure |2] shows more rapid variability due to the in- 
creased sampling with respect to the weekly-averaged LAT 
light curve. Overlapping time intervals of the 18-month 
light curves in 7-rays and in the R band were used to cal- 
culate the discrete cross correla tion function (DCCF), fol- 
lowing lEdelson & KrolikI (Il988h . In Figure |9] the DCCF is 
plotted with error bars estimated by a Monte Carlo method, 
takin g measurement errors and data sampling into account 
dPeterson et al.l ll998). The correlation strength and lag were 
computed by fitting a gaussian profile to the DCCF between 

*^ |http ://users.utu. f i/kani/lm/| 

"^^ [http : / /james . as . arizona . edu / ~p5mith/Fermi| 
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Fig. 9. — Discrete cross correlation function (DCCF) between tlie 7-ray 
and optical light curves of Fig. |2] 

lag -20 to -1-20 days. The result was a correlation intensity of 
0.17 ± 0.09 (corrected for the effect of measurement noise) 
and a lag of —5 ± 5 days, where negative lag means 7-rays 
lagging the R band. The choice of lag range for the fit, 
from -15, -1-15 to -30, -1-30, as well as bin size, 3 or 4 days, 
only had a small effect on the estimated time lag of the peak 
and on the correlation strength, with the latter ranging from 
0.14 to 0.17. The rather weak long-term correlation between 
gamma-rays and optical tends to disfavor a one-zone syn- 
chrotron plus synchrotron self-Compton (SSC) model, since 
in this model one would expect electrons with approximately 
the same energies to make both the optical and 7-ray emis- 
sion, and might favor an external radiation Compton (ERC) 
or multi-zone synchrotron /SSC model (e.g.. lAharonian et al.l 
I2009t iBottcher et al.l 12009). The optical-UV variability ob- 
served during the 48 days of the multiwavelength campaign 
is not correlated with the weaker X-ray variability seen by 
Swift -XRT and RXTE nor with the non-variable 7-ray emis- 
sion observed in the same period by the LAT (Figure|4|i. The 
lack of correlation might support the scenario where both syn- 
chrotron and IC photons contribute to the X-ray emission. 
Variability produced by the highest energy electrons would 
be smeared out by the IC component, and the mix of pho- 
tons from the low and high energy component would dilute 
the correlation with the synchrotron, or synchrotron plus ther- 
mal, optical-UV emission. This would explain why variability 
in the optical/UV is greater than in X-rays. Moreover if the 
synchrotron and IC emission components come from phys- 
ically distinct regions, they would also weaken the correla- 
tion. No evidence for a contribution by accretion disk radia- 
tion during this low activity state of the campaign was found, 
although this might contribute to the lack of clear correlation 
between the X-ray and optical-UV emission. The recent dis- 
covery of a luminous Ha emission line in the op tical s pectra 
of BL Lac (.Vermeulen et al. 1995: Corbett et al.ll 19961 boOOt 
ICapetti et al.ll2010l) . as likelv a product of the broad fine re- 
gion photoionized by the disk radiation, can suggest a more 
relevant role played by the disk during major outburst with a 
high 7-ray Compton dominance (as seen in the EGRET 1997 
outburst). 

The averaged SED from the MW campaign(2008 Aug. 20 
- Sep. 9) is reported in Figure[TO] For the first time the multi- 
frequency emission for the eponymous blazar source BL Lac 
is mapped from radio to GeV 7-rays during a low-activity 
state. The relevant archival observations of BL Lac are also 
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Fig. 10. — 2008, Aug. 20 - Sep. 09 simultaneous SED of BL Lac assembled with data of the coordinated campaign, and some archival data from the literature 
for comparison. The data are modeled with a single zone pure-SSC (top panel), two-zone pure SSC description (center panel), and a single zone SSC plus ERC 
descriptions (bottom panel). In the two-zone pure SSC description the blue/dark dotted line represents the compact, faster and high-energy emission component 
while the blue/dark dashed-dotted line is the larger and diluted zone emission component. In the SSC plus ERC model the blue/dark is the ERC emission (disk 
emission negligible here) component and the blue/dark dashed-dotted line is the s ynchrotron and SSC emission component. Past published and ai'chival data are 
taken from Bottcher & Bloom l2000): .Ravasio et al.n2002l) :iB6ttcher et al.' ^200?): 'Albert et al. (2007); Raiter i et al.i 12009, ). and fit models for two of these past 
epochs (the 1997 outburst campaign and the 2005 campaign, Ravasio et al. 2002; Albert et al. 2007) are represented as orange/clear long-dashed lines in each 
panel. 



reported in this figure for comparison in addition to our data. 
The difference in 7-ray integrated flux with respect to the out- 
burst state of 1997 measured by EGRET is about 15 times 
lower 

The XRT spectral results are in agreement with the 



RXTE results in the common energy range 2-10 keV, and 
with the usual SED properties of ISP blazars. In fact the con- 
firmation of a flat X-ray SED ( vF^ ~ const) of moderate 
flux level with both instruments in such range, makes us more 
confident the X-rays are from a transition region between the 
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synchrotron and inverse Com pton energy components. In ad- 
diction, as reported in Section lTTl the whole 0.3-lOkeV range 
spectrum accumulated during the campaign by Swift -XRT 
showed a broken power-law best fit, i.e. a steeper soft X- 
ray spectral slope which should be the tail of the synchrotron 
emission and a flatter slope above 3 keV, which should be the 
first part of the IC energy component. 

The LAT 7-ray spectrum accumulated in 18-months is re- 
ported as well for comparison (open/orange circles represent- 
ing non simultaneous data in Fig. [TOl l. 

A single-zone synchrotron self-Compton (SSC) descrip- 
tion, usually successful for TeV blazars, has been used to 
model also the SEP of thi s eponymous BL Lac in the past 
dRavasio et al J 120021 l2003h . and is thus the first step in the 
attempt to reprod uce its broadband emission. This model 
(iFinke et al.l2008l) is applied to the SED of BL Lac (blue/dark 
thick line on the top panel of Figure [Toli. The energetic blob 
is modeled with Doppler factor V ~ T — 2Q, comoving ra- 
dius R = 5.3 X 10^ cm, (which is consistent with a 2 hour 
variability time observed in the source on previous occasions), 
and tangled magnetic field of intensity B = 0.33 G. The elec- 
tron distribution has the form of a broken power law, so that 

A^e(7) CX -f-P^ fov-frrnn = 600 < 7 < -/break = 1-4 X 10^, 

and A^(7) oc '^^P" for -jbreak < 7 < Imax = 10'' where 
7 — E/{mc^) is the electron Lorentz factor The best fit 
electron indices were found to be pi = 2.0 and p2 = 3.8, 
which is consistent with electrons in the fast cooling regime. 
It is interesting to note that 7„,i„ is well constrained, for a 
given a blob size, by the X-ray observations, and thus it is 
impossible for this model to reproduce the radio points. This 
means that, if this model is correct, the emitting region must 
be closer to the black hole than the radio photosphere. In 
this model the Poynting flux power and the kinetic power 

in relativistic electrons are = 1.5 x 10'''^ erg s^^ and 

JL.Z <^ 

Ljgj = 2.2 X 10"''^ erg s^^, respectively, implying an equipar- 

tition fraction eg ^ U'^/K = ^^tV^et = 0-068. This is 
well below the equipartition value, which may be argued as a 
reason to disfavor this model. In this case the Doppler factor is 
considerab ly higher than the radio-variability Doppler factor 
reported in Hovatta et al.l (2009), although they evaluated that 
from the estimated timescale of each radio flare component 
and the known brightness temperature, referred therefore to 
a much larger region beyond the radio photosphere than this 
model suggests. 

Since the one-zone SSC model has difficulty reproducing 
the variability of this source and giving reasonable energet- 
ics (eB substantially below 1), as further step we attempt a fit 
with a two-zone SSC model (blue/dark thick line in the cen- 
tral panel of Figure [Tol i. using two emitting blobs via SSC 
emission: a compact and faster emission region and a larger, 
slower and diluted region accounting also for the radio-band 
emission. The energetic blob responsible for most of the emit- 
ted power and originating in a limited part of the jet is mod- 
eled with Vi 10, Ri = 9.0 X 10^^ cm, Bi = 0.45 G. 
The radio-band and hard X-ray emission is reproduced us- 
ing a much larger emitting region characterized by parameter 
values V2 = 6.5, R2 = 1.6 x 10^^ cm, and B2 = 0.02 
G. Doppler factors of the two zones are found to be closer 
to the avera ge value (V = 7.3) found in previous years 
by iHovatta et al.l (12009) than the one-zone SSC model pre- 
sented above. The kinetic partial differential equation of this 
model describes the evolution of the particle energy distribu- 



tion A^(7, t) after the injection of shock-accelerated electrons, 
with a rate equal to a power law with an exponential cutoff: 
(5(7) oc 7-^6"^/'*'"'°^, [cm~^ s~-^], where 7 = E/{mc^), 
between ^min and ^ mpx- More details on the numerical 
model can be found in (ICiprinill2008l 120 1 Ol) . For the compact 
region (1) we obtain 7^\^„ = 3.0 x 10^ and -/'ikL = 7.0 x 10^ 
with soft injection index p*^^^ = 3.5, and L^.^^jj^^ — 6.1 x 10^^ 
erg s^^, ijgjfi-) = 2.2 x 10*^ erg s^^, implying an equiparti- 
tion fraction ~ 2.8. Considering the larger diluted region 
(2) we obtain 7,*,J]„ = 600 and 7i,?L = 3.0 x, lO'' with in- 
jection index ^ 2.4, and i^ftja) ^ ^-^ x lO'^^ erg s"\ 
^jet(2) = 7.9 X lO**^ erg s~^, implying an equipartition frac- 

tion e^' = 0.020. In this case for the two regions the to- 
tal €b = 2.16, a value slightly closer to energy equipartition 
than the one-zone model. The hypothesis of two synchrotron 
components present at different distances from the nucleus of 
BL Lac and accounting for the optical and th e X-ray spec- 
tra respectively, was suggested in the past in iRavasio et alj 
(2003). In our case X-rays could be considered produced 
from both the emitting regions (central panel of Figure [TOb , 
and this can explain the lack of evident correlation between 
UVOT and XRT/RXTE light curves. On the other hand the 
fact that the variations in 7-rays are not seen, or not resolved 
by the LAT, contrary to the optical-UV variable flux, tell us 
that the two-zone pure SSC model does not explain the multi- 
frequency light curve of Figure |4] observed in BL Lac during 
the campaign, although the energetics is closer to (and above) 
equipartition with respect to the previous single-zone model. 

SSC emission with the relevant addition of external ra- 
diation Compton (ERC) emission was used to explain the 
large Compton dominance of BL La c during the 1997 
7-ray outburst recorded by EGRET (Bloom etal. 119971: 
Madeiskiet"anil 999: Bottcher & Bloom 2000; Ravas io etal J 
20021: iBottcher & Reimer 2004|). An accretion disk emit- 
ting component, and emission lines were also observed sev- 
eral times in this blazar jV ermeulen et al. 19951: ICorbett et all 
[19961 12OOOI: ICapettietalfcolO ). The SSC plus ERC hybrid 
model, where X-ray emission is produced via the SSC mech- 
anism a nd the GeV 7-r ay emission is produced by ERC, was 
used bv lMadeiski et al.l (11999) for this source. This motivates 
its use to see if it is able to reproduce also the 2008 low state 
of BL Lac. A model fit to this SED is shown as the blue/dark 
thick line of the bottom panel of Figure [10] For this model 
fit an equilibrium version of the time-dependent jet model re- 
ported in Bottcher & Chiang (2002) was used. A power-law 
distribution of electrons, (3(7) = Qo7~^> with lower and up- 
per cut-off 7mi„ = 1.1 X 10'^ and 7m a a: = 2.0 X 10^, respec- 
tively, and injection index p = 2.85 is injected into a compact 
region of radius i? 3 x 10^^ cm, moving relativistically 
along the jet, oriented at an angle 6'obs = 3.8° with respect to 
the line of sight, with bulk Doppler factor 2? = 15 and tan- 
gled magnetic field intensity B = 2.5 G. The instantaneous 
electron injection is self-consistently balanced with particle 
escape on a time scale <osc = VcscR/c, where t^csc = 60 is 
a free parameter, and radiative cooling through synchrotron, 
SSC and ERC emission. For the ERC emission, both the di- 
rect accretion disk radiation field and accretion disk emission 
reprocessed in the BLR are taken into account, although for 
our fit, the BLR-reprocessed radiation strongly dominates the 
radiation which is directly from the disk. We assume an ac- 
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cretion disk luminosity of Ldisk = 6 x 10"'^ erg s^^, and an 
emitting region distance from the black hole of — 0.1 pc. 
The BLR is represented as a spherical shell of reprocessing 
material with radial Thomson depth tblr — 0.01, located at 
a distance tblr = 0.21 pc from the central black hole. Al- 
though the model fits the radio data well, the assumed distance 
from the black hole within the BLR implies that the blazar re- 
gion is well within the radio-band photosphere. The numeri- 
cal model self-consistently evaluates the energy content in the 
resulting equilibrium electron distribution, and compares this 
value to the magnetic-field energy density. This model gives 
L^.f^) = 4.75 X lO'^^ gj.g g-i ^(e) ^ 3 2 x lO^a, and ratio 
cb = 1.48. This is closer to equipartition than the one- and 
two-zone pure SSC models described above. Considering the 
equipartition ratio es the differences between the two-zones 
SSC and the ERC scenario do not appear particularly signif- 
icant because of the uncertainties in the estimate of this pa- 
rameter, but the uncorrected variability can still be explained 
more easily in the ERC scenario. 

The parameters used for the three SED models that are 
shown in Figure [TO] are in agreement with the limit on the 
magnetic field obtained from the VLBA observations, while 
the Doppler factors used in the two-zone SSC are in better 
agreement with the value extracted from radio-band structure 
historical observations. The magnetic field and Doppler fac- 
tor could be different for the 7-ray emit ting region than for 
the radio core region described in section [l!2l since the 7-ray 
emittion region could be from a smaller region than the ra- 
dio core, closer to the black hole. A ERCh-SSC model can 
produce more complex SEDs and less correlated variability 
patterns, partly because of variability in the external radiation 
field. The emission region could move through an inhomo- 
geneous radiation field, when it leaves the BLR. Synchrotron 
and SSC depend on the magnetic field, while the ERC com- 
ponent does not. Therefore, magnetic field fluctuations may 
also produce uncorrected variability patterns between opti- 
cal - X-rays (SSC) on the one hand and 7-rays (ERC) on the 
other. 

4. CONCLUSIONS 

This is the first broad-band multifrequency campaign in- 
vestigating the low activity state of BL Lac. This eponymous 
source is also pecuhar, sometimes showing a rather complex 
SEP profile, a behavior sim ilar to an FSRQ ( Vermeulen et alj 
11995^'Mad ejski et alJll999 ). and being detected also by VHE 
Cherenkov telescopes hke MAGIC above 200 GeV. We ob- 
served a 7-ray flux (E > 100 MeV) during the 48-day mul- 
tifrequency campaign (planned intensive campaign, PIC) of 
2008 about 15 times lower than the level measured by EGRET 
during the 1997 outburst, and lower than the 1995 EGRET 
observations. This value and the weekly light curve of this 
campaign interval indicates a quiescent 7-ray state. The aver- 
age LAT spectrum accumulated during the first 1 8 months of 
survey showed evidence for curvature. The single photon of 
maximum energy detected by the LAT during these months 
was measured at 70 GeV, with a spectrum that could extend 
into the VHE range. The weekly light curve is characterized 
by at least two mild flares on the 1-week timescale (in the 
middle of 2009 April and the end of 2010 January). 

During the PIC the low luminosity and non-variable 7-ray 
state did not correspond to the lowest luminosity state for 
near-lR and optical emission. No evidence for clear corre- 
lation between the low-energy emission component (radio to 



UV bands) and the high-energy emission component (X-ray 
to 7-ray bands) was found. Since one would expect correlated 
variability between these components in a one-zone SSC sce- 
nario, a fit with a two-zone SSC model was attempted. This 
model as well is not able to explain the variable optical and 
UV emission seen by Swift corresponding to the non-variable 
7-ray emission during the coordinated campaign. Therefore a 
hybrid SSCh-ERC model was better able to reproduce the SED 
being closer to energy equipartition and providing a possible 
explanation of the uncorrected variability. 

This result tells us that the SSCh-ERC model is valid for 
BL Lac also during low brightness states (strengthening the 
quasar-like, FSRQ, view of this blazar). The two-zone SSC 
model might still be a possibility for this blazar (as sug- 
gested, for example in iRavasio et alj 120031) . Other blazars 
like S5 0716- H71 appear to have evidence for tw o synchrotron 
components ("Giommi et al."2008|: IChen et al .''2008). In 3C 
454.3, the prototypical FSRQ, a two-zone synchrotron sce- 
nario is inferred by the varia bility of the d ust IR emission 
bump resolved by Spitzer tOgle et alj|2010l) . Therefore, as 
suggested in the past, BL Lac, a rather special object with a 
possible Seyfert-like nucleus, might be considered a miniature 
and real FSRQ (Vermeulen et al. 1995; C orbett et al. 199^ 
iMadejski et aTir999l:ICorbett et al. 2000: lCapetti et alJI2 010). 

A multi-frequency VLBA observation performed during 
this campaign demonstrated that the bright, innermost part of 
the radio-mm jet ("the core") cannot be described by a single, 
uniform, self-absorbed, synchrotron emitting blob, but rather 
an inhomogeneous region with spatially changing spectral 
properties. The synchrotron self-absorption turnover is de- 
tected in the VLBA data which allows us to constrain the av- 
erage magnetic field in the core region to be less than 3 G, as- 
suming a particular Doppler factor from previous radio single- 
dish observations. SED model values are in agreement with 
this limit, although this is not strictly necessary since the mod- 
els assume emission from a region closer to the black hole 
than the radio core. 

Another interesting discovery was that during the cam- 
paign, BL Lac continued to show a broken, concave X-ray 
spectrum, as seen by the Swift -XRT This was not confirmed 
by RXTE , probably because of the differing energy range 
and worse sensitivity of this instrument, but the photon in- 
dex Tx ^ 2 is in agreement with this view. We note that 
the RXTE spectrum was mostly consistent with the XRT one. 
Both instruments observed an X-ray spectrum with small day- 
to-day X-ray photon index variations during this low-activity 
state. The concave spectrum and lack of spectral index vari- 
ation is believed to be a typical signature of intermediate en- 
ergy peaked blazars, and we find it exists even during non- 
flaring states. 
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